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The frustrated classical antiferromagnetic Heisenberg model with Dzyaloshinskii-Moriya (DM) 
interactions on the triangular lattice is studied under a magnetic field by means of semiclassical 
calculations and large-scale Monte Carlo simulations. We show that even a small DM interaction 
induces the formation of an Antiferromagnetic Skyrmion crystal (AF-SkX) state. Unlike what is 
observed in ferromagnetic materials, we show that the AF-SkX state consists of three interpene¬ 
trating Skyrmion crystals (one by sublattice), and most importantly, the AF-SkX state seems to 
survive in the limit of zero temperature. To characterize the phase diagram we compute the average 
of the topological order parameter which can be associated to the number of topological charges 
or Skyrmions. As the magnetic field increases this parameter presents a clear jump, indicating a 
discontinuous transition from a spiral phase into the AF-SkX phase, where multiple Bragg peaks 
coexist in the spin structure factor. For higher fields, a second (probably continuous) transition 
occurs into a featureless paramagnetic phase. 

I. INTRODUCTION 

Twisted modulated magnetic textures due to an 
antisymmetric spin exchange interaction, termed the 
Dzyaloshinskii-Moriya interaction, have attracted much 
interest mainly after the experimental observation of 
non trivial magnetic configurations, called magnetic 
Skyrmion lattices, which have important potential tech¬ 
nological applications^. Recently, a strong evidence of 
the formation of a Skyrmion crystal (SkX) state was ob¬ 
served in metallic ferromagnet MnSi^. This state consists 
of a triangular lattice arrangement of Skyrmions and can 
be visualized as a superposition of three non-equivalent 
spirals states (each characterized by one wavevector 
In systems with a square lattice structure, the 
SkX phase arises from the competition of the ferromag¬ 
netic and the Dzyaloshinskii-Moriya (DM) interactions 
and it is stabilized by a magnetic field and thermal fluc¬ 
tuations. Recently, Okubo et.al have shown that it is 
also possible to stabilize the SkX phase in the isotropic 
Heisenberg model in the triangular lattice with strong 
nearest-neighbor (ferromagnetic) and weak next nearest- 
neighbor (antiferromagnetic) interactions^. 

In this paper, based in analytical approximations and 
Monte Carlo simulations, we show that in a pure antifer¬ 
romagnetic frustrated system it is possible to stabilize a 
phase closely related to the SkX phase described above. 

The novelty resides in the fact that the state found in the 
present paper consists of three interpenetrating Skyrmion 
lattices (one by sublattice) and, most importantly, this 
antiferromagnetic SkX (AF-SkX) state survives in the 
limit T ^ 0. Unlike what happens in the ferromagnetic 
case where it is necessary to include large values of D/| J|, 
in the antiferromagnetic case small values oi D j J can 
stabilize the AF-SkX state. To identify and characterize 
this phase we compute a chiral order parameter which de¬ 


termines the density of Skyrmions present in the lattice. 
Due to the discrete nature of the model, in the region 
of large magnetic fields a series of discrete jumps in the 
order parameter is observed. 

The paper is organized as follows: In Sec. [iTlwe present 
the Hamiltonian and study the classical ground state so¬ 
lutions at zero magnetic field and zero temperature based 
on the spherical approximation. This analysis shows that 
the ground state has a three-fold degeneracy which could 
be a source of exotic ordered states that might be real¬ 
ized under applied fields, e.g. various types of states 
where multiple wave vectors coexist (so-called multiple-q 
states). Sec. IHII contains our Monte Carlo results. We 
found a quite rich low temperature behavior of the sys¬ 
tem as the magnetic field is varied. The system goes from 
a spiral low field phase to an antiferromagnetic Skyrmion 
lattice phase at higher field. Finally, increasing the mag¬ 
netic field further, the system enters into a featureless 
paramagnetic phase. We pay particular attention to the 
undoubtedly most interesting phase which is the antifer¬ 
romagnetic Skyrmion lattice phase. It has the particu¬ 
larity of being composed of three interpenetrating and 
shifted Skyrmion lattices that realize on each sublattice 
of the triangular lattice. The MC study is complemented 
in Sec. [IVlby a phenomenological analysis and in Sec. [Vl 
we analyze the stability of the Skyrmion phase and briefly 
discuss the relation of the present study with the square 
lattice analog. We conclude in Sec. [^with a summary 
and discussion of our results. 


II. MODEL AND SPHERICAL 
APPROXIMATION 

Let us consider the antiferromagnetic Heisenberg 
model on the triangular lattice in a magnetic field, with 
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the Hamiltonian given by 

n = Sr-Sr,+D6r-{Sr X Sr')-h'^S^ (1) 

(r,r'> r 

where the spin variables are unimodular classical vectors, 
|Sr| = 1, J > 0 and D are respectively the antiferro¬ 
magnetic and DM couplings, 6t = (r' — r)/|r' — r| is a 
unitary vector pointing along the axis, rp the 

sum over nearest-neighbors (NN) couplings on the tri¬ 
angular lattice with primitive translation vectors of the 
direct lattice ei = (1,0) and 62 = (1/2, \/3/2) and h the 
strength of the magnetic field along the 2 ; axis. The DM 
interaction is chosen in the direction of the (5r vectors, so 
as to give rise to spiral spin states with the spins lying in 
a plane perpendicular to the propagation vector—. 

The magnetic phase diagram for the model defined 
by Eq.([T]) with D = 0 has been discussed iniid^. At 
T = 0 and for zero magnetic field, the system orders 
in a 120° three-sublattice magnetic structure described 
by the wave vector k = (47r/3,0). In a magnetic field 
the classical energy is minimized for spin configurations 
constrained by the magnetization of each triangular pla- 
quette: 


Sa = h/(3J) . (2) 

This constraint leaves undetermined the orientation of 
the spin plane and sublattice directions inside that plane. 
This degeneracy persists up to the saturation field hg = 
9J. However, due the order-from-disorder effect^ this 
degeneracy is lifted by thermal fluctuations, which se¬ 
lect a kind of states (collinear and coplanar) over the 
non-coplanar ones at low temperature. The resulting 
phases depend on the magnetic field strength: for low 
field a coplanar so-called Y-state, with one spin pinned 
in the negative z direction and two canted up, at exactly 
h/J = 3 the pseudo-plateau state (at M = 1/3) with the 
collinear configuration up-up-down and in the high field 
region a coplanar canted version of it, which smoothly 
interpolates with the fully polarized state at h/ J = 9. 

Now we turn-on the DM interaction with the orien¬ 
tation along the nearest neighbors bonds (Fig. [1]). As 
a first step we investigate the magnetic order of the 
Hamiltonian Eq. m at zero magnetic field and zero 
temperature by means of the spherical approximation^. 
Within this scheme, instead of imposing the local length 
constraint |Sr| = 1, one imposes a milder condition, 
= N where N is the number of lattice sites. 
With this softer constraint, the model Hamiltonian o 
can be diagonalized by a simple Fourier transformation 
^r,a — Here, the index a = z is the 

spin component, a = 1,2,3 is the sublattice label and 
r and k denote the position and pseudo-momentum re¬ 
spectively. The Hamiltonian becomes 



FIG. 1. (Color online) (a): Triangular lattice. The green 
arrows indicates the primitive translation vectors of the di¬ 
rect lattice ei = (1,0)? ^2 = (l/2,v^/2). Red arrows the 
Dzyaloshinskii-Moriya vectors. The sites with labels , ri ...r 4 
indicate the sites involved in the calculation of the local chi¬ 
rality. (b) Area of the triangle to compute for the discretized 
Skyrmion number (Eq. Q). 
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where the 9x9 matrix M(k) 


M(k) = 


0 mi2 mis 
mi2 0 ^32 

^13 ms2 0 


( 4 ) 


and the 3x3 matrices mah depend on J, D and k with 
explicit expressions given in Appendix 
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FIG. 2. (Color online) Minimum eigenvalue ie™^’^(k) in the 
spherical approximation for D/J = 1/2. For any value of 
D/J > 0 we find that the ground state corresponds to a 
triple-^ state by sublattice. 


k 


( 3 ) 


Since we are at T = 0, the ground-state is obtained 
from the minimum eigenvalue of the matrix 
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M(k) (Eq. dH)). In particular, we find that for any 
value of DjJ > 0 there are three minima (see Fig. 13 ) 
suggesting the presence of multiple-q states. In this con¬ 
figuration the ordering wavevectors k* appear along the 
directions of the nearest-neighbor bonds with module 


ik*i 


D V3D2 ^3 
7 4 J2 “ 4J3 


32 


+ 0 ( 



In particular we are interested in the non trivial triple- 
q state, which consists of a superposition of three spirals 
(each characterized by one wavevector k*), that corre¬ 
sponds to the Skyrmion lattice phase. Within the spher¬ 
ical approximation, performed at both zero temperature 
and zero magnetic field, we find it conceivable that, at 
very low temperatures, one may neglect all wave vectors 
other than the three critical modes k* and the uniform 
k = 0 mode. Henceforth, we expect that it maybe pos¬ 
sible to find ordered states characterized by the number 
of wave vectors. In order to elucidate if multiple-q states 
could be stabilized at finite temperature and magnetic 
field we explore the behavior of the system by means 
Monte Carlo simulations in the next Section. 


III. MONTE-CARLO SIMULATIONS 


Monte Carlo simulations are performed based on the 
standard heatbath method combined with the over¬ 
relaxation method. Periodic boundary conditions were 
implemented for N = L‘^ site clusters with L = 12 — 60 
and L = 72,96 for some calculations. The system is 
gradually cooled down from high temperature. A run 
at each magnetic field or temperature contains typically 
0.1 — 1 • 10^ Monte Carlo steps (MCS’s) for initial relax¬ 
ation and twice MCS’s during the calculation of mean 
values. 

The different phases can be easily identified from the 
observation of real-space spin textures (e.g. see Fig. [31). 
We have also computed the average of various physi¬ 
cal quantities in order to identify precisely the different 
phases and the corresponding transitions between them. 
We first calculate the standard magnetization, magnetic 
susceptibility and specific heat. 


M = 


C = 



NT‘^ 



(6) 

(7) 


We also introduce the discretized Skyrmion number XQ 
and the total chiralit y xl^ defined as 


XQ = ^ (8) 

(f) 


where = ||(Sr„ - SrJ x (Sr, - SrJ||/2 is the lo- 

cal area of the surface spanned by three spins on ev¬ 
ery elementary triangle Vi^Va^Tb (see Fig. □o). Here 
= Sr^.(Sr„ X SrJ is the so-called local chirality 
and r^, ri ^ are the sites involved in the calculation of 
Xq (see Fig. [T]). It is important to note that for slowly 
varying spin configurations XQ Xl would coincide. In 
the case at hand the use of XQ turns out to be more ef¬ 
fective because in the discrete case it approximates much 
better the directed area of the sphere surface, therefore 
providing a much better measure of the Skyrmion num¬ 
ber. We have computed the discretized Skyrmion num¬ 
ber by sublattice and so the sites Vi^Va^rb run in one 
sublattice. However, the value of the local chirality with¬ 
out sublattice distinction is also an important quantity 
within the context of the anomalous Hall effect. Indeed, 
electrons moving in the antiferromagnetic background of 
our system, and within the framework of the adiabatic 
approximation, would feel an effective flux per plaquette 
which is given by xl—^^ computed between sites of the 
three different sublattices. 

With the help of XQ we can detect the Skyrmion phase 
as a function of the magnetic field as shown in Fig. |4| 
(Center). At low magnetic field h < yg = 0 due 

to the spiral configuration. Then for h*(T) < h < 
a non zero valor of xq indicates the presence of a wide 
AF-SkX phase. For instance, in this region and for a 
system size of A = 48^ the discretized Skyrmion num¬ 
ber in each sublattice Xq/^c ^ 12/Ac = 1.5 x 10“^ 
(Ac = A/3) which is in perfect agreement with the num¬ 
ber of Skyrmions that can be directly observed from the 
MC snapshots (in Fig. [3] only a portion of the total lat¬ 
tice is shown). We also plot xl together with XQ showing 
that Xq is a better measure of the Skyrmion number. We 
observe that the discretized Skyrmion number is almost 
constant in the complete region and decay in few steps 
as a consequence of the discrete nature of the Skyrmion 
patterns. In this region of the phase diagram all the 
spins have ^(-component greater than zero and its xy- 
components look like an array of 2D-vortices, a pattern 
that can be interpreted as the remnant of the Skyrmion 
lattice. The stepwise decrease of the Skyrmion number 
should be associated with the finite size of the lattice, 
which, due to the periodic boundary conditions imposed, 
can accommodate only definite numbers of Skyrmions. 
In the continuum limit we expect that this number will 
smoothly go to zero. In order to test this idea, we have 
performed an approximate (continuum limit) computa¬ 
tion of the Skyrmion radius as a function of the magnetic 
field, which leads to a smooth decrease of the Skyrmion 
density (see Sec. IB. 

In all the MC simulations the spin configurations have 
a three sublattice structure, so in order to understand 
better the underlying structure, we plot the patterns 
by sublattice. A, B and C. As an example, in figure [3] 
we show representative spin configurations for L = 48 
(only a small region of the entire lattice is shown to illus¬ 
trate the spin texture) at three different magnetic fields: 
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(i) Spiral 


(ii) AF-SkX 


(iii) high-field (iv)AF-SkX phase: three sublattice 










FIG. 3. (Color online) Snapshots by sublattice for L = 48, D/J = 1/2, T/J = 9 x 10 ^ and h/J = l{i),2A{ii) and 6.2(m). 
The green box indicates which sublattice is ploted (A, B or C). In figure (iv) we show all three sublattice together. 


hlJ= (i)l, {ii)2A and (iii)6.8. 

We have also calculated the static spin structure factor 
in the reciprocal lattice to identify the Bragg peaks that 
characterize the different spin-textures. The perpendicu¬ 
lar and parallel (to z) components S±{q) and 5'||(q) are 
defined as 

5±(q) = ^(1 E + I E (10) 

r r 

r 

where () means the averaged MC configurations. In fig¬ 
ure [5] we show the intensity of the spin structure factor 
for the patterns found (figure [3]): spiral or single-q (Top) 
and AF-SkX or triple-q (Center) and high field phase 
(Bottom). 

(i) Spiral phase (FiglfiH): the spin structure consists 
of three interpenetrating spirals on each sublattice, 
A, B and C (see figure[T]). Each one is characterized 
by one of the three possible ordering wave vectors 
k* obtained by the Spherical Approximation (Eq. 
[5] and figure [2]) . Both the xy and components are 
characterized by the same wave vector as can be 
seen from the figure O 

(ii) Antiferromagnetic Skyrmion lattice phase (Eigifib): 
here the stable phase corresponds to three SkX 
phases, one on each sublattice. Each of the SkX 
phases is a superposition of three spirals character¬ 
ized, both the xy and z components, by the three 
wave vectors k*. The complete superposition of 
these three sublattice structures is shown in Eiglfiji. 

(iii) High field phase (Eiglfib) which in some areas shows 
a vortex-like structure (Eigifib). There, the xy com¬ 
ponents of the spins form vortices while the z com¬ 
ponent is always positive. This phase is estab¬ 
lished near the right edge (higher field sector) where 


the chirality order parameter, both xg, xl show 
smaller plateaux 


The main result of our study is summarized in the 
phase diagram shown in Eig. [4] (Bottom), which was ob¬ 
tained for D/J = 1/2 (other values of this ratio lead 
to a similar phase diagram). The AE-SkX phase is sur¬ 
rounded by the Spiral and paramagnetic phase and the 
lines separating the phases were obtained in different 
ways, locating the peaks in the specific heat (figure [6]), 
the peaks in the magnetic susceptibility (Eig. [4|(top)) 
and using XQ the order parameter (figure [4|(center)). 

A very interesting behavior is observed when the sys¬ 
tem goes from the AE-SkX phase to a region where the 
XY components of the spins form a vortex-like struc¬ 
ture. In all the AE-SkX region, as we increase the mag¬ 
netic field, the Skyrmion number does not change. To 
increase the magnetization, the Zeeman term tends to 
align the spins, then in a big region of the phase dia¬ 
gram, antiferromagnetic exchange and DM interactions 
favor the AE-SkX phase with a specific topological num¬ 
ber (the number of Skyrmions). If we continue increasing 
the magnetic field, the system (being discreet and in a 
phase with a given number of Skyrmions) cannot increase 
the magnetization without paying a high energetic cost, 
thus what is observed is that the AE-SkX phase changes 
the number of Skyrmions. This result is shown in figure 
[4|(Center) where both the discretized Skyrmion density 
and total chirality density decrease by steps. 

In the next Section we use a phenomenological- 
variational analysis to compute the dependence of the 
Skyrmion radius with the magnetic field. 
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FIG. 4. (Color online) Top: Magnetization M (blue) and 
magnetic susceptibility vs. magnetic field h for L = 48, 
D/J = 0.5 and T/J = 9 x 10“^. Center: Discretized 
Skyrmion number XQ and total chirality xl, both quantities 
per sublattice, vs. h/J forL = 48, D/J — 1/2. Red triangles, 
green rhombus and purple triangles indicate the discretized 
Skyrmion number A, B and C. Blue dots show the total chi¬ 
rality. Bottom: Complete h — T phase diagram. The limit of 
the obtained by the peaks in the specihc heat and from the 
changes in the chirality order parameter. 


IV. PHENOMENOLOGICAL ANALYSIS OF 
THE SKYRMION LATTICE SPACING 


Although in the MC data the number of Skyrmions for 
a given system size, and hence the lattice spacing varies 
with the applied magnetic field by discrete jumps, we 
suspect this behavior to be a finite size effect. In the 
thermodynamic limit one can expect a lattice spacing 
continuously varying as a product of the competition be¬ 
tween the applied magnetic field which tends to increase 
the magnetization by making bigger Skyrmions and the 
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FIG. 5. (Color online) Intensity plot of the static spin struc¬ 
ture factor for L = 96 for the spiral phase (Top); AF-SkX 
phase (Center) and disordered phase (Bottom). 



FIG. 6. (Color online) Specihc heat vs. temperature for 
h/J = 1.4, 3 and 6. The positions of the peaks are used 
to locate the lines separating the different phases. 


DM interaction for which a hxed Skyrmion size is op¬ 
timal. To model this competition, we follow the same 
idea developed by Jung Hoon Han et. in the case of 
a ferromagnetic system. Of course here we have a three 
sublattice Skyrmion pattern so the phenomenological de- 
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scription of a free energy in terms of spin orientations 
that we present has to be interpreted as representative of 
an average free energy for one of the three sublattices. 

Consider the local spin orientation ( 6 >, 0) of a single 
Skyrmion depending on the local coordinate (r, 0) as 0 = 
if — Tij2 and 0 = 0{r). The total free energy of a single 
Skyrmion then reads 


Fsk = 47 r J J rdr 

sin^( 6 >) 


IdO \ 2 ^ 

2 *+“) 


+ - 


4^2 


— r]{cos{0) — 1 ) 


( 12 ) 


where 77 = h/(2J) and k = D/(2J). The function 0(y) 
minimizing the functional ([T 2 ]) describes the magnetiza¬ 
tion distribution in an isolated Skyrmion. To understand 
the observed transitions phenomenologically, we calcu¬ 
late the energy of a isolated Skyrmion using the following 
trial function 


0{r) = TT, r < Ro 

1 


= TT- 


r—Rp 
R—Rp 


1 


^ R-Rp 


Rq < r < R 


= 0, R < r 


(13) 

(14) 

(15) 



R 


FIG. 7. (Color online) Dependence of the free energy density 
given by Eq. (fT6|) as a function of the Skyrmion radius R for 
two values of the Skyrmion core: Ro = 0.5,1 = 1) 


V. SKYRMION PERSISTENCE AND THE 
IMPORTANCE OF THE LATTICE 


In the definition of 6 >(r), the parameter Rq > 0 is the 
minimum radius simulating the central spin in the con¬ 
tinuous limit, R > Ro is the radius of the Skyrmions 
and 6 is a parameter to be computed self-consistently. 
The parameters R, b are fixed by the condition of mini¬ 
mum, {dFskx/dR^dFskx/db} = {0,0}. The Skyrmion 
lattice observed in our simulations (in each sublattice) 
correspond to a close-packing of individual Skyrmions in 
a triangular lattice. We then write the total free energy 
of the Skyrmion lattice state as 

where L is the size system of the system. In Fig. Owe 
show a few plots of Fskx as a function of parameter R, 
i.e. as a function of the Skyrmion radius. We observe 
that the behavior of the optimal Skyrmion spacing as a 
function of the magnetic field varies very slowly for small 
values of the magnetic field and then behaves much more 
abruptly when approaching the high field paramagnetic 
phase. In a fixed finite size sample, where the Skyrmion 
lattice has to be commensurate with the sample size, this 
translates precisely in a large zone of the phase diagram 
in which the Skyrmion number is fixed and then start to 
produce jumps that appear more and more abruptly. 


It is well known that Skyrmions are topological defects 
with an infinite life-time in the continuum and with a 
very long persistence time in ferromagnetic lattice mod¬ 
els. Here we are going to argue that despite having an 
apparently more fragile three sublattice structure, the 
Skyrmion persistence is also quite noticeable. To do this 
we start from a configuration, like figure [3] in the AF-SkX 
phase, obtained by annealing at fixed magnetic field. In 
Fig. [HI we show an example of the MC time evolution 
(the curves are the average of 500 copies) of such initial 
state in our simulations after a quench to zero magnetic 
field. Note that unlike what is observed with the nor¬ 
malized magnetization M (x0.05 to use the same scale 
for all quantities), the total chirality xl and discretized 
Skyrmion number by sublattices XQ exhibit greater stiff¬ 
ness extending over several decades of time. 

To complete the analysis we have also performed MC 
simulations in the non-frustrated square lattice where 
in contrast to the case of the triangular lattice, = 
(r' — r)/|r' — r| is an unitary vector pointing along the 
axis ei = (1,0) and 62 = (0,1). Our simulations do not 
show any evidence of an antiferromagnetic Skyrmion lat¬ 
tice in the square lattice case. One possible explanation 
for the absence of this phase is geometrical frustration. 
At the microscopic level, on the square lattice, spins have 
the tendency to form a two sublattice structure, instead 
of three for the triangular lattice. On the other hand. 
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MC steps 


MC steps 


FIG. 8. (Color online) The Monte Carlo time evolution 
of the magnetization, discretized Skyrmion number XQ ^-^d 
total chirality xl dehned by Eqs. calculated at 

T = 9 X 10“^ J for sizes of L = 24,36,48, 60. The initial 
state corresponds to the AF-SkX phase (Fig. [3|). Each curve 
corresponds to an average of 500 independent time evolutions. 


stable for a wide range of magnetic fields and couplings. 

Observing the evolution of the topological order pa¬ 
rameter as the magnetic field increases, which we de¬ 
fine as the chirality per sublattice and which can be in 
turn associated to the topological charge or number of 
Skyrmions, we observe a clear jump, indicating a discon¬ 
tinuous transition from a spiral phase into the AF-SkX 
phase, where multiple Bragg peaks coexist in the spin 
structure factor. For higher fields, a second (probably 
continuous) transition occurs into a featureless paramag¬ 
netic phase. 

We also observe that after a magnetic field quench, 
the Skyrmion number persists for longer MC times as 
one should have expected, making apparent the robust¬ 
ness of the AF-SkX state. The relation of this AF-SkX 
state with the 3-sublattice structure of the lattice studied 
in the present paper is contrasted with the square lattice 
case, where such a state does not show up. Finally, we 
have also shown that doped electrons that would evolve 
in this background, would feel a fictitious non-zero mag¬ 
netic fiux, as depicted in the central figure of[4l Inter¬ 
estingly, the total chirality, which is directly related to 
the local flux felt by the doped electrons, also remains 
quite constant in agreement with the fact that the AF- 
SkX structure is remarkably stable for a wide range of 
magnetic fields, and this of course can have very inter¬ 
esting consequences in the context of the anomalous Hall 
effect. 
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Appendix A: Matrix Elements for Spherical 
approximation 

The explicit expressions for the 3x3 matrices of the 
matrix 0 are given by 

TOl2(k) = [Jf2'^],mi3(k) = [jf3^],TO32(k) = 

where 

Jif = (l + 7tk + 7l.kM“'' 

2 

+K ^(ei - 7i* 1,62 - 72,ke3)^e“'^^ (A2) 

7=1 

< = (i+7r,k72,k+7r.k)<^“'' 

2 

+K ^(62 + 7lk72.ke3 - 7i,kei)^e"'^^ (A3) 

7 = 1 

'^^2 = (1 + 7l,k72,k + 72, 
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+K ^(^3 + 7i,k72.ke2 - 72,kei)^e“'^^ (A4) 

7=1 

where 7 = (1,0), 62 = (1/2, y^3)/2), 63 = ei-e 2 and 


7 i k = and 72 , 1 , = e-**"'( 2 e 7 -e 2 )_ 
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